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Abstract. For the description of the transport of electrons across a quantum dot, 
which is tunnel coupled to leads at different chemical potentials, it is usual to assume 
that the total Hamiltonian of the composite system of the leads and the quantum dot 
is the sum of three contributions: That of the leads (nonintcracting electrons), that 
of the quantum dot and a third one, the "tunneling Hamiltonian" , which reflects the 
possibility that electrons can move from the leads to the quantum dot or vice versa. 
The text aims at a mathematically clear derivation of such a separation. I will start the 
discussion with the total Hamiltonian of the system acting on a many-electron wave 
function, including the attractive interaction between nuclei and electrons as well as the 
repulsive Coulomb-interaction between different electrons. Indeed, a natural separation 
of the total Hamiltonian in the described form will be obtained. An analysis of the 
tunneling Hamiltonian shows that the electron-electron interaction yields contributions 
to it which represent the correlated tunneling of two electrons at the same time. For 
the derivations it was useful to introduce a map called "antisymmetric product" . In an 
appendix I show possible exact representations of the total Hamiltonian (with arbitrary 
potential V(r)) obtained by the use of the antisymmetric product. 
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Tunneling Hamiltonian 



2 



1. Hamilton operator 

The Hamilton operator of a crystal acts on an A-electron wave function 

iff: D N -> C, (1) 

where D = M 3 x {|, \.} , D N = D x D x . . . x D. It contains the kinetic energy, the 
attractive force between the nuclei and the electrons and the repulsion between different 
electrons. It can be written as 

h 2 N 2 

^ = + £ V ^ + £ uTh^- ( 2 ) 

i=l i<j ' 1 3 1 

The position of the nuclei is considered as constant in time within this equation pp. I 
summarized in "V(r)" the potentials of the bare nuclei located at the lattice points of 
the crystal. A wave function as in ([ID can be viewed as a set of 2 N maps 

VV,...,^ : {M. 3 ) N -XC, 

(n, • • -,r N ) V'«i,...,^( r i» ■ • • ' r ^)- 

The Hamiltonian leaves the spin- variables cr, G {t, 1} unchanged and acts on each of 
the 2^ components. 

The above shape of the Hamilton operator is quite general. It should be valid even 
if the nuclei are not located at the points of a regular lattice and also, for example, if one 
considers two finite crystals separated by vacuum. The changes enter into the potential 
V(r). 

Lateron in the text it will turn out to be convenient to have an extra notation for 
the contribution of the electron-electron repulsion to the total Hamiltonian. I define 
thus 



i<j 



2. Antisymmetric wave functions 

2.1. Slater determinants 

A wave function of N fermions is required to be antisymmetric. Slater determinants 
[1] provide a way to construct antisymmetric wave functions of N electrons by the 
use of one-electron wave functions: For a system of one-electron wave functions 
•0i : D — y C, i — 1, . . . , N, one can define 

as 

ipi® ...®ip N (x 1 ,...,x N ) := — = ^ sS^WVVM^i) • • -^{n){xn)- 

I used "xj" to denote elements of D, " SV' is the set of all permutations of the set 
{!,..., N}. 
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(For the set of permutations of any finite set the following rules hold: Any 
permutation can be represented as a product of subsequent "transpositions" . These are 
those special permutations which exchange only two elements and act, apart from this, 
as identity. The representation of an arbitrary permutation as product of transpositions 
is not unique, but the property "even" or "odd" of the number of needed transpositions is. 
Therefore, one can assign a sign to the permutation, sgn(n) = ±1, in case the number 
is even/odd. This sign is multiplicative, sgn(TCTr') = sgn(7r)sgn(ir'). The "Leibniz 
formula" uses permutations to express the determinant. ) 

At any rate, the Slater determinant is an antisymmetric wave function. For any 
two (quadratically integrable) wave functions ifj,ip' : D N — > C one defines the scalar 
product pQ as 

<V#'> := E • • • E / C...,^!, • • • ^N^'^ajn, . . . ,T N ). 

2.2. Normalization condition 

The norm is defined by || ip || 2 = (iplip). For simplicity, the prefactor l/v^VT was 
included in the Slater determinant. It ensures that, if the functions ipi,---,' l pN form 
an orthonormal system, then the norm of their Slater determinant is one. To fill 
the normalization condition with life, I tried to clarify its relation to the particle 
density. In the case of N distinguishable particles represented by a wave function 
ip(ri, . . . , tat, 0tv) one could speak about a probability of finding particle 1 with spin 
0i within a volume V\ C M 3 and . . . and particle N with spin 07V within a volume 
Viv C M 3 . The probability might be defined as 

1 2 



dn... dr N \ip{rx,ax,...,r N ,a N )[ 
'Vi Jv N 

Correspondingly, one would define the number of particles which can be found in a 
volume V C K 3 as 



El/ +/ + 2 / }lVW 2 (n,r 2 )| 5 

CT 2 KJVxV c JV c xV JVxVJ 



I assumed the special case N = 2 and used the notation V c = R 3 \ V for this. One can 
then represent 



+ / +2 = + 

>VxV c JV c xV JVxV JVxR 3 JR 3 xV 

The way of defining the probability and representing it in an alternative form can be 
generalized to the case of arbitrary N. If the one-electron wave functions ipi, . . . , ipN 
contributing to a Slater determinant form an orthonormal system, then one obtains 
indeed the particle density 

N 

p( r ) = EE i^( r ' a )i 2 

i=l a 

which one would expect intuitively. 
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2.3. Antisymmetric product 

A Slater determinant can be viewed as an antisymmetric product of wave functions 
ipi : D — > C. A generalization is possible: For 

if):D m ^C, <p:D n ^C 

one can define 

if) ® <p : D rn+n -> C 

by 



^ <8) tp(x u . . . ,x m+n ) := a - — — V sgn(-K M ) 

\ (m + n ! 

V v y M:Mc{l,...,m+n},\M\=m 

^(.•E-itmWi ' • ' ' ^im(»i) )V- ; '(^'7rM(m+l) ' • • • ' ■^^M{m+n) ) • 

I used "|M|" to denote the number of elements of M. The permutation 7T/vf is the one 
which counts at first the elements of M in increasing order and then the elements of its 
complement M c :={l,...,m + n}\ M. More precisely: 

7r M (l) = xi, . . . , 7Tm(w) = x m , Tc M (m + 1) = y 1 , . . . , 7Tm(w + n) = y n 

if M = {xi, . . . , x m }, M c = {yi,...,y n } in increasing order. By the property 
sgn(7iM)sgn(7iM c ) = (— \) mn one can show the relation 

ip®p = (-i) m7 V <g> -0. 

If if) and <y? are antisymmetric, then also ip ® (p. This can be seen by applying the 
"antisymmetrizer" [2] to if) and ip. 

The product is bilinear and associative: 

if) ® (<£ ® A) = (V» ® <p) ® A. 

A symmetric representation is given by 

[V> ® (v ® A)] (xi, . . . , x m+ „+i) = 



m!n!Z! v-^ \ 

/ — ; rrrr sgn{7r M ,N,L) 

[m + n + /)! 

v y M,N,L: MUNilL={l,...,m+n+l},\M\=m,\N\=n,\L\=l 

V ; ( a; 7rA/ i ]v,L(l)' ■ ■ ■ J X ^M,N,L{m)) 
( f{ X 'TM,N,L(m+l)j * * * J S 7TM,iV,i(m+n)) 
A( a '7r Mi] v,L(m+n+l) > • • • j " C 7i-m,]V,l(™+™+0 ) ' 

where 7Tm,n,l counts at first the elements of M, then the elements of A/", and finally the 
elements of L in increasing order. 

As a consequence, expressions like u tf>i <$...<%> i/)n" (without any brackets) are 
well-defined for ipi : D mi — >■ C. Indeed, if all Wj = 1, then the product is the Slater 
determinant defined above. 
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3. Vector space 

3.1. Vector space of the leads 

About the vector space of the many-electron states of the leads I assume that it 
is generated by single-electron wave functions via Slater determinants and linear 
combinations of these. This is in agreement with the concept that the many-electron 
states in a large crystal can be specified by saying which "one-electron levels" are 
occupied. The Fermi-Dirac distribution was derived in Ref. [1] by the use of this idea. 
According to the Born-von Karman boundary condition one has two such one-electron 
wave functions for every "allowed" wave vector A; in a primitive cell of the reciprocal 
lattice. I denote the one-electron wave functions by 

i)i ka :M 3 x {|,|} 

The index "Z" denotes the lead, ipik{r) is a Bloch wave function depending only on 
a space- variable. A further band index might be included but is not essential to the 
purpose of this work. The number of one-electron levels is proportional to the size of 
the crystal pQ. The ipik a form an orthonormal system in the quadratically integrable 
functions D — > C. For every natural number n <\R\ one can define 

V n (R) := lin{ipi ® . . . ® ip n : ipi, • • • ,ip n € R} • 

(With u lin(S)" I denote the set of all linear combinations of elements in S.) V n (R) 
is a linear subspace of the set of all antisymmetric and quadratically integrable maps 
D n — > C. The subsets M C R with \M\ = n form an orthonormal basis of V n (R) via an 
identification 

M — ipx <g) . . . <g) ip n , where ipi, . . . , ip n 

is a list of the elements of M in some previously chosen and fixed order. The set of linear 
combinations of Slater determinants is really larger than the set of Slater determinants. 
In general, the solutions to the eigenvalue equation Hip = Eip are linear combinations 
of an infinite number of Slater determinants. Any antisymmetric function can be 
approximated by a sequence of linear combinations of Slater determinants [5], while 
there is in general no way to approximate a linear combination of Slater determinants 
by a sequence of Slater determinants. The ground state energies of two-electron atoms 
have been calculated numerically by the use of the Hamilton operator, Eq. (T5]) [5 J . The 
good quantitative agreement with experiments indicates that this rather intransparent 
Hamiltonian is indeed the right one as long as relativistic effects are ignored. Moreover, 
the Hartree-Fock equations have been derived from this Hamiltonian. They contain an 
exchange term which has been made responsible e.g. for screening effects. 
Essentially in analogy to Ref. [3] I define the vector space V(R) as 

V{R) := V {R) x V 1 (R) x . . . x V W {R), 
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where V (R) := C. Since the V n (R) are Hilbert spaces via their scalar products 
(.|.), V(R) is also a Hilbert space. V n (R) can be viewed as a subset of V(R) by the 
identification 

^eV n (R) corresponds to (0, . . . , 0, ip, 0, . . . , 0) G V{R). 

The element (?/>o, • • • , i>\R\) £ ^(-R) can be written as + • • • +^|-R| • The scalar product 
is 

l«l 

((tpO, ■ ■ ■ ,^\R\)\(fO, ■ ■ ■ , f\R\)) = 

Moreover, one has the map 

<g> : V(R) x -> 

i+j<\R\ 

The elements of V"(i?) can be viewed as linear combinations of subsets of R, in the same 
way as the elements of V n (R) can be viewed as such linear combinations. The element 
1 G C = Vo(R) represents the empty set. With this interpretation one can write for 
M, N C R: 

M®N = ±M UN in case M n N = 
and zero otherwise. 



3.2. Approximate Hamiltonian 

The Hamiltonian ll H Rl1 of the reservoirs is given by Eq. (T5]) in case one replaces the 
potential V(r) appearing there by the potential of the reservoirs "V (r)". (V(r) = 
^ R (r) + V®(r) with V^(r) the potential which enters the Hamiltonian of the quantum 
dot.) There is little hope that H R is indeed an endomorphism of the Hilbert space 
V(R), i.e., in general the image H R (V(R)) is not a subset of V(R). However, I assume 
that for states ip G V(R) with "sensible" electron numbers H R tp is very close to being 
an element of V(R) and that in this sense V(R) is a good choice of the vector space on 
which the Hamiltonian operates. By the use of the orthogonal projection 

Pv(R) : lin (V(R) U H R (V(R))) -)• V(R), 

which can be defined by the condition 

i>-pv(R^ 1- V(R), 

one obtains the approximate Hamiltonian pv( R )H R . This is indeed a hermitian operator 
V(R) — > V(R) and defines the mathematical model of the reservoirs. 
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3. 3. Vector space of the composite system 

The construction of a vector space V(M) can be realized for any finite set M of one- 
electron wave functions. In case M is an orthonormal system the subsets of M form an 
orthonormal basis of V(M). I assume that the vector space of the many-electron states 
of the isolated quantum dot is given by V(Q) where Q is a finite orthonormal system 
of one-electron wave functions. Analogously to the above arguments, I assume that for 
the "most relevant" states if of the quantum dot H®(p is very close to being an element 
in V(Q). Hence, a good approximate Hamiltonian might be pv(Q)H® . 

I feel forced to assume that R U Q is still an orthonormal system. This assumption 
is in agreement with previous tunneling theories [7] . Going even beyond this, I demand 
that 

H R (V(R)) C VjQ 1 ), 

H Q (V(Q)) c VW). (4) 

U R ± " is the orthogonal complement of R in the square-integrable functions D — > C. 
V(R ) is the set of all linear combinations of Slater determinants of elements in R ± , 
in the same way as V(R) is the set of all linear combinations of Slater determinants 
of elements in R. Finally, Vi^R^) is the topological closure of V^R^), i.e., the set of 
all infinite and convergent linear combinations of Slater determinants of elements in 
R- 1 . The demand is that, in case one applies the Hamiltonian of the reservoirs to any 
wave function in the chosen vector space of the leads, then the result is a wave function 
for whose representation as an (infinite) linear combination of Slater determinants of 
one-electron wave functions exclusively elements in Q 1 - are needed, and vice versa. This 
should be approximately fulfilled if the distance between the sub-systems is large. 

(Remark: V(R) and VlR- 1 ) are not orthogonal since C is contained in both of 
them. Still, they are somehow antisymmetric complements since for all ip £ V(R),<p £ 
V(R ) :|| ip <g> ip || 2 =|| if> \\ 2 \\ f \\ 2 - One might take this equation as a definition 
of an antisymmetric complement "a.c. " and think about the question whether indeed 
V(R) a - c - = VXRPj. ) 

As vector space of the many-electron states of the composite system of reservoirs 
and quantum dot V(R U Q) should serve. In case of large distance between the two 
systems the assumption makes sense, for sure: The eigenstates of the composite system 
can be expected to be just the product of the eigenstates of the separate systems. 

The Hamiltonian, Eq. (J2J), acts on V(R U Q) like 

H(4>oi if>i, ■ ■ ■ , i)\mjQ\) = (0, Hifjx, . . . , Hip lRuQl ), 

i.e., in particular iJ0 = 0. As the final model Hamiltonian I take 

Pv(ruq)H: V(RUQ)^V(RUQ). 

The vector space V(R U Q) is a realization of the "tensor product" of the 
vector spaces V(R) and V(Q), which is defined in mathematics in a purely formal 
way. The antisymmetric product is a realization of the corresponding bilinear map 
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"Cg> : V(R) x V(Q) —> V(R U Q) u . The purely formal tensor product has been used in 
Ref. j2] to construct the vector space of a composite system of two sub-systems. 

4. Additive separation of the Hamiltonian 

For a start, I want to consider the action of the total Hamiltonian on an antisymmetric 
product ip ® ip G V(R U Q) with ip G V m (R), <p G V n (Q), 1 < m < \R\, 1 < n < |Q|. 
For this I write the electrostatic potential V(r) caused by the nuclei entering the total 
Hamiltonian, Eq. (]2]), as V(r) = V®{r)+V R {r) where V s (r) is the electrostatic potential 
entering the Hamiltonian of the corresponding sub-system. The product ip ® ip can be 
written as 



1>® ( P = \ 7 r~\T sgn(ir M )S WM [ip * <p] , 

\ [m + n)l z — ' 

V v ' Mc{l,...,m+n}:\M\=m 

where ip * ip : D m+n — > C is defined by 

ip * (p{x u x m+n ) = ip{x x , x m )(p(x m+ x, x m+n ) 
and where for any permutation 7r G Si, I G N, and any function 

a : D l C 
the function S n a : D l — > C is given by 

S n a(xx, ...,xi)= a(x^, . . .,x n Q)). 

This representation is convenient since the Hamiltonian commutes with the operators 
0, 



H[ip®<p] = J- — — V sgn(7r M )SK M H [ip *<p] 

y (m + n ! z — ' 

V v ' M<z{l,-,m+n}:\M\=m 

With a little bookkeeping [TJ one can write 

# [ip ® <£>] = (tf'V) ®<p + ip® {H Q <p) 



/ 771 71' x — 

+ \/7 T~\T 2^ sgn(-K M )S 7rM T mtn [ip * , (5) 

1/ (m + n)\ r *— ' , , , 

where T mjn is defined by 

rp ._ T-i(Q) I T-i(fl) I rp(int) 

J-mjTi • -"-771,71 ^ -'■m,?! ^ - 1 m,n ' 

T S } H r l, O-l, • • • , r m+n, Vm+n)] ■ = 
m m+n 9 

— - — - e z 



It-. _ r .\ f a ( r l'°"l' • • ■ > r m+n,<7m+n)> 

«=1 j=m+l 



T S H r l> fl> • • • , r m+n, 0"m+n)] - 
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j &1 ) ■ ■ ■ ) '"m+ra ) C"m+n J j 



i=l 



and where Tm)i has the corresponding definition. 

For the above alternative representation of H [if) <g) </?] the assumption m, n > 1 was 
made. The equality, however, is correct also in the cases that m or n or both of them 
are zero. (For example, ifx&ip = if)ip in case if) or ip is only a number; TJ^o = 0.) Roughly 
speaking, the terms in Tm*n reflect the repulsion between electrons in the contacts on 
the one hand and electrons in the quantum dot on the other hand. The terms in T^l 
represent the attractive interaction of the electrons in the reservoirs on the one hand 
and the nuclei in the quantum dot on the other hand. 

Since pv(ruq)H was chosen as the model Hamiltonian, the projection pv(ruq) should 
be applied to both sides of Eq. (jSj). By the demanded complementarity PJ, one obtains 
the equation 

p V (RUQ)H [ip ® If] = {j) V{R) H R i)) <g> ip + V> ® (p V (Q)H Q if) 



/ 7TI Tl x. — *\ 

+Pv(ruq) \ ( m + n )\ 2^ sgn(TT M )S VM T min [if) * tp] . 

V ^ Mc{l,...,m+n}:\M\=m 

Proof: 

I used that Pv(ruq) [{H R ip) ® <p] = (p v ^ R) H R ip) ® (p for any if) G V m (R),<p G 
V n (Q), m, n > 0. Thus, I ought to show that 

(H R if) — pv( R )H R if)) <g> <p JL V(RUQ). 

According to the complementarity-condition H R if) is a (maybe infinite) linear 
combination of Slater determinants of functions in Q L . Since R is a finite subset of 
Q- 1 , one can assume that all of the Slater determinants appearing in the expansion have 
the shape 

lf)l ® ■ ■ ■ lf> io <g> lf> io +l <8> • • • ® V'm, 

where ipi E R if i < i , ipi-LR otherwise. If one applies the projection py(#) to such 
a Slater determinant, then one obtains zero in case iq < m. The Slater determinant 
remains unchanged in case io = m. Thus, H R if> — pv(R)H R if) is a linear combination of 
such determinants with i < m. For any ip G V n (Q) and any such Slater determinant, 
however, one gets 

if)x® ...if) io ®if) io+1 ® ...®ip m ®<p 1 V(RUQ). □ 

For any two endomorphisms A G End(V(R)),B G End(V(Q)) (End(V) is the 
set of linear maps V — > V) there is a unique well-defined endomorphism "A £g> B Ll 
G End(V(R U Q)) which has the property: 

A <g> B(ip ® ip) = (Aip) <g> for all if; G p G V(Q). 
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The proof of this statement is straightforward and perfectly analogous to the proof of 
the corresponding statement about formal tensor products. One can define: 



H t ot 
Ho 



PV(RUQ)H, 

(p v(R) H R ) g> l v(Qh 

lv(fl) ® (PV(Q)H Q ), 



and finally 

-Ht := H tot — Hr — Hq. (6) 
All of these are hermitian operators on V(R U Q). 

5. Creation- and annihilation operators 

By the use of creation- and annihilation operators [3] H? can be expressed in a compact 
form. For any normalized element ip G V\(R U Q) one can define the creation operator 
of -0 as the map 

4 : V(RUQ) ^V(RUQ), 

ip i— > ip (g) </?. 
The annihilation operator of the one-electron level ^ is defined as the adjoint operator 
of the creation operator, := ^c^j . 

6. Expression for in terms of creation- and annihilation operators 

So far, a representation of the total Hamiltonian 

Htot = H R + Hq + H T 

was obtained, where can be called the contribution of the leads, Hq is the 
contribution of the quantum dot and is the tunneling Hamiltonian. The main 
part of this text aims at a description of Ht- Ht leaves the particle number unchanged, 
since the other contributions as well as the total Hamiltonian are doing so; Hq% = 0. 

For the description of the linear operator H T one can use any basis of V(R U Q). 
Independently of whether Hr and Hq are diagonal in the Slater determinants given by 
the subsets of R U Q one can thus use these Slater determinants for the description. 
For example, if different elements in R } if yi,...,y n are n different 

elements in Q and if z±, . . . , z m+n e RUQ are m + n different elements, then one obtains 
according to Eq. (jSJ) 

(zi <g> . . . ® z m+n \H T [xi <g) . . . (2) x m g) yx <8> . . . <g> y n ]) = 

- y y y y 

(m + n)\ ^ ^ ^ ^ 

v ' Mc{l,...,m+n}:\M\=m Tr€S m+n rGS m r]£S n 

sgn(n M )sgn(n)sgn(T)sgn(r]) 

( z n(l) * ••• * ^7r(m+n) | 

Sn M Tm,n [ x r(l) * ■ ■ ■ * ^r(m) * Vr)(l) * ■ ■ ■ * 2/r?(n)] )• (7) 



Tunneling Hamiltonian 11 
6.1. Contribution of the electron- electron interaction to 

First, the contribution of to the matrix element is considered. Hence, I assume 

that m, n > 1. For any 1 < i < m, 1 < j < n, the permutation G S m+n is defined by 
(myself and) the conditions: 

• Kij(l) = i, Kij (2) = j. 

• Kij(3), . . . , Kij(m + n) is an increasing list of the remaining numbers. 

For any sequence of one-electron wave functions one can write 

e 2 

-Vi * ... * v m+n = 

In-rji 

S Kij { [H {mt) (v Kij{1) * v Kij{2) )] * v Kij{3) * ... * v Kij(m+n) } . 

For this, the operator of the electron-electron interaction defined in Eq. ((3]) was used. 
Moreover, I define now 

w x := xi, . . . , w m := x m , w m+1 := y 1 , . . . , w m+n := y n , (8) 

moreover the permutation " r ■ r\ " G S m+n by 

• ( r • n)(l) := t(1),1 < m. 

• (r • r]){l) := r]{l — m) + m, I > m + 1. 

Inserting all of this into Eq. ([7]), the contribution 

Mc{l,...,m+n}:|M|=m 7rg5 m+n rG5 m r?eS n i<mj>m+l 

sgn(ir M )sgn(7r)sgn(T)sgn(r]) 

( z mr M (T-ri)- 1 (T-r})Ki j (l) * • • • * ^ttttmCt"-')) -1 (T-v) K ij (m+n) \ 

H {mt) [w( T . v ) KtjW * tW( T .„)« ii (2)] * W( r .^) Klj ( 3 ) ... * W( r .^) Klj ( m+n) ) (9) 

of the electron-electron interaction to the matrix element is obtained. The fact that 
the operators SV' conserve the scalar product, (S^iplS^ip) = (ip\<f), and the relation 
n ii were used for this. (In particular: S w/ = S^r-i.) 
For all 1 < i < m, m + 1 < j < m + n and f , g G {1, . . . , m + n} , f ^ g , I 
define the set Sf°J o ° as the set of all tuples 

(M,7T,T,T], 

with the property: 

(r-r])Kij(l) =i , 
(r-T])Kij(2) =j , 
7rvr M (r • rj)' 1 ^) = fo, 
TrK M {r ■ 77) _1 0'o) = go- 
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The set of all tuples (M,tt,t,T), over which the sum goes is the disjoint union of 



the sets Sf°f . The sum (19]) turns into 



fogo 

J. 11C OUlll ipil LUlllO 111LU 

sgn(irir M {r ■ i])- 1 ) 



E E 

OoJo,/o,So) (M,7r,r,r),(iJ))esf 



90 
0J0 



(m + n)! 



JJ (^ M (r.^)-l(/)kz)- (10) 

There is at most one permutation u n{°j°" £ 5 m+ „ with the property that the 



product 



WO 

«^*0 JO 



is nonzero and with the property that 

The value of the product is one in this case, //there is a permutation 7r^°J o ° with the 
demanded properties, then the inner sum of expression ( flOl) turns into 



sgn{4:f Q ) (z fo *z go \H^[ Wio *w jo )), 



«0J0 

since the number of tuples (M, ir, r, 77, (z, j)) in ^{°j° with the property tttt m (t ■ rf)^ 1 = 
ir{°!-° is then (m + n)\. Swapping the position of fo and g has the effect that the sign 
switches, so one obtains for f < g the equality 

(M,7r,T, J? ,(ij))GS /oso US 9o/o 
V WO WO 

= ^9<4To) ( z fo ® l# M) K ® «>*]> (12) 
in case there is a permutation tt^ with the demanded properties (fTTT) . 
The sign can be represented as 

s 9 n (^iojo) = <*1 ® ■ ■ ■ ® ^m+nlc^cLc^o^io [™1 ® • • • ® ™m+n])- 



The equality f fT2|) can be rewritten as 



e sg " ( ^; ) :' r ' ) (**.^->k.» j j> 

(M,7r,r,r;,(ij))GS /o90 US 9o/o 
V WO WO 

(^TrTTMCr-r?)- 1 ©!^) 

^OjO 

(*i <g> . . . <8> z m +n|(2/ ® %,l# (int) No ® U> io ]) 

C L C l C ^o C <% kl ® ■ ■ ■ ® Wm+n})- 
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Written in this way, the equality is correct for all 1 < i < m,m + 1 < j < m + n 
and all 1 < fo < go < m + n, no matter whether a permutation with the properties ffTTl) 
exists. The sum (El) f=(TT0l)') turns into 



m m+n 

{ Zl ® ...®z m+n \Y^ E E (zf^ZglH^lwi^Wj]) 

i=l j=m+l l<f<g<m+n 
[w 1 <g) . . . ®W m+n }). 

Going back to Eq. ([7]), one can write the contribution of the electron-electron interaction 
to the matrix element 

(zi <g> . . . <g> z m+n \H T [x\ ® . . . <2) x m ® y x ® . . . ® y n \) 

as 

m n 

( Zl ® . . . ® z m+n \J2J2 E ( z f®*g\ H{ira) [xi®Vi]) 

i=l j=l l</<<?<«i+n 

= (zi ® . . . <8> z m+n \H^ nt) [xx ® . . . <g) x m <g) yi ® . . . ® y n ]), 
where I used the definition 

H T Ut) --= E E (*®P'\H {mt) a®P) clc\,c p c a . (13) 

a,a'eR/3,/3'eQ 

is a hermitian operator and is the contribution of the electron-electron interaction 
to the total tunneling Hamiltonian. 
Remark: 

jjiyn-t) alternative representation 

4 mt) = l/2EE E (7®5|# {mi W/3> c\c%c a . 

To verify this, one has to show that in case 7 G R and 5 G R the matrix element 

(j®6\H<Ma®P) 
vanishes. One writes this as 

( [(H {mt) -H R ) + H R ] 1 ®5\a®(3) 



and uses that if fi (7 ® 5) G V^Q -1 ) according to the complementarity-condition 01]). 
Because /3 G Q, the contribution of vanishes. A close look and another application 
of the complementarity yield that also the contribution of H^ nt ^ — H R vanishes. 

6.2. Contribution of the potentials to Ht 

Next, the contribution of Tm-n to the matrix element of Eq. (JTj) is determined. For this 
I assume m > 1 and write in the same way as in the previous subsection (Eq. [8]): 

X T (i) * ... * X T ( m ) * * ... * y v (n) = W(r. v )(l) * ... * W( T .^)(i). 
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In analogy to the previous subsection, one can write for any 1 < % < m and any sequence 
v i, . . . , v m+n of one-electron wave functions 

V Q (ri) [vi* v m+n ) = S Ki [(V Q (r)v Ki ^)) * v Kt{1) * ... * v Ki{m+n) ] , 

where the permutation Ki e S m+n , which is defined by the conditions 

• Ki(l) = i, 

• Ki(2), . . . , Ki(m + n) is an increasing (arbitrary) list of the remaining numbers ^ i, 

was used. One obtains the contribution 

sgn(im M (T ■ n) -1 ) 



E E 



(m + n)! 

M,ir,T,r] i=l v ' 



(^7r7r M (r-r,)- 1 (r^) Kl (l)|^ /Q (^)^(r- }7 ) Kl (l)) 
m+n 



of the potential V® to the matrix element (j7|). For every 1 < i < m and every 



l</o<' m + ^I define the set E[° as the set of all tuples 



'0 



(M, 7T, T, T), l) 

with the property that 

(r ■ r])Ki(l) = i , 
7T7r M «i(l) = fo- 

The set of all tuples over which the sum goes is the disjoint union of the sets E{°. Thus, 
the sum can be rewritten as 

E v sgn(7T7r M (T ■ n)' 1 ) 
^ (m + n)! 

»0,/o (Af,7r,r,n,j)GSf o ° 

(z f0 \VQ(r)w i0 ) 

X\( Z TnT M {T-q)-^(l)\wi)- 

For fixed Zo, fo, there is at most one permutation u tt{°" with the properties 

• IWv°(oK> ^° 

• *£(io) = /o. 

If there is such a permutation, then the number of all tuples (M, tt, r, n, i) G E^° with 
the property that 

7171 M (r -7])- 1 = 7r{° 

is (m + n) ! . The inner sum has in this case the value 

sgn{4° ) (z h \V^{r)w l0 ) 
= (zi <g> . . . <g> z m+n \ (zf \V Q (r)w io )cl fo C WiQ (wx ® . . . ® u> m+n )). 
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The representation of the inner sum by the second line is correct even in the case 
that there is no permutation with the properties of Thus, the contribution of the 
potential V® to the matrix element of Eq. (J7|), 

(zi <g> . . . ® z m+n \H T [xi® . . .®x m ®yi® . . .® y n }), 

is 

(zi <g> . . . <g> z m+n I 

m m+n 

i=i /=i 

[si®...®a; m ®yi®...®y„]) = 
(^i <£)...(£) z m + n \ 

E E (^i^coa) 4 c « 

aGi? /SGfflJQ 

[xi g> . . . g> x m g> yi <g> . . . <g> y n ]). 

The sum of the contributions of the potentials V® and V R to the tunneling 
Hamiltonian is 

+ E EE ^H*w> 4c,. (M) 

se{R,Q} aes pes 

I used the letter U S U to denote the two sub-systems; is the the complementary 
sub-system of S; V(r) is the sum of the two potentials, V = V R + V Q . 
Remark: 

has the alternative representation 

4 V) = E E E MV 8 (r)a) 4c a . 
Se{R,Q] aes peRuQ 

To see the identity, one can use the assumed complementarity (@J and write for 
a e R, (3 e Q: 

(a\V R /3) = ((H R -cA)a\/3) = (-cAar|/9) 

= (a\-cAp) = (a\-{H Q -V Q )(5) 

= (P\v Q ay, 

where I used the abbreviation c := — 

2m 

6.3. Expression for Ht in terms of creation- and annihilation operators: Summary 

The tunneling Hamiltonian, defined by Eq. 02]), is the sum of the contributions of the 
electron-electron interaction, Eq. ffT3]) . and that of the potentials, Eq. f JT4|) : 

H T = E E («'®/3'l^ (mt) «®/3) 44, C £ Cq 
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+ E E H^VK) el*, 

Se{R,Q} a,a'eS 

+ E EE \(«\nr)P) ofa. 

S£{R,Q} a£S /3eS 

6-4- Example 

I would like to illustrate a possible effect of the additional contributions , Eq. (fl3|) . 
to the tunneling Hamiltonian. If I assume that a, a' 6 R, /3, fl' G Q with 

a (r,cr) = V/fc(r)5j,(cr), 
a'(r,a) = if>i'k>(r)6<i(cr), 
/3(r,a) = ip t (r)5 t (a), 
(3'(r,a) = Vl (r)S i (a), 
then the corresponding contribution to the tunneling Hamiltonian reads 
.2 f [ ,,rn( r )*i ) m(r)(pt{r')i>i>k^(r , y jj 



-e J J drdr ^— — c^c^c^. (15) 

The term annihilates an electron with spin "down" in the leads as well as an electron 
on the quantum dot with spin "up" . At the same time, it creates an electron with spin 
"down" on the quantum dot and an electron with spin "up" in the leads. The term 
represents a spin-flip process and simultaneous scattering of electrons in the leads. The 
coefficient which is obtained in case I' = I can be expected to be larger compared to the 
case I' I (I is the lead index). 



7. Outlook 



By taking the reduced density matrix and applying the projection operator technique 
[2], it might be possible to approach the transport problem with two leads at different 
chemical potentials even when terms as in the example of the last subsection are 
included. A complete (diagrammatic) analysis of the kernels as in Ref. [10] would 
be desirable, but not necessary for perturbation theory of lowest order in the tunneling 
Hamiltonian. The terms in Ht of the form 

(a\V(r)P) c^cp, aeS,(3eS, 

give rise to energy conserving single electron tunneling seen by perturbation theory of 
lowest order in the tunnel coupling, e.g. Ref. [10]. Analogously, one might expect that 
from the terms of the form of the example (1151) one obtains transition rates 

iV(t-4) 

whose value is sensitive to the question whether there are electrons and holes with 
appropriate spins in the leads, such that the spin-flip- and scattering process obeys the 
energy conservation. This might have interesting consequences for the current across 
the quantum dot, especially if / ^ V . 
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8. Conclusions 

I considered the Hamiltonian of an electronic system which consists of two weakly 
interacting sub-systems. The starting point was the many-electron Hamiltonian 
obtained from the Born-Oppenheimer approximation. Moreover, it was assumed that 
the wave functions of the states of the composite system can be constructed by the use 
of the wave functions of the states of the sub-systems in a natural way (antisymmetric 
product). The assumption of weak interaction became manifest in a "complementarity 
condition" 01]). It turns out that the total Hamiltonian has a natural additive separation 
into three contributions: For each of the sub-systems one obtains a contribution which 
is essentially given by the Hamiltonian of the corresponding sub-system. 

The contribution which expresses the interaction is called tunneling Hamiltonian. 
An analysis of it showed that it can be represented in a natural way in terms of 
annihilation- and creation operators. The contributions to the tunneling Hamiltonian 
are classified according to their origin: One kind of terms is caused by the electron- 
electron interaction, while another kind of terms is due to the interaction of the positive 
nuclei and the electrons. 

The considerations are very general: The inner electronic structure of the sub- 
systems was not relevant. Even if the electron-electron interaction can be taken into 
account in the sub-systems by an effective noninteracting Hamiltonian, it is still expected 
that in the tunneling Hamiltonian the terms caused by the electron-electron interaction 
appear. 

In my mind and in the text I used the picture that the systems are distant 
and separated by vacuum. As an example, the conditions for the derivations are 
approximately fulfilled also in the case of impurities in a metal. The sub-systems are 
then the impurity and the metal. The overlap of the one-electron wave functions can be 
expected to be small if the level generated by the impurity is localized. Apart from the 
transport across a quantum dot, one possible way to see whether the spin-flip scattering 
processes contained in the tunneling Hamiltonian yield sensible results would be to apply 
them to the problem of a resistivity minimum in metals doped with magnetic impurities 
[I], explained first in Ref. [8] by spin-flip scattering processes. The formal similarity 
of the problem with the transport across quantum dots led to the prediction of the 
zero bias resonance known as Kondo resonance [111 El [9] . Spin-flip scattering processes 
with even a similar formulation as in the example ( 1T51) have been taken into account 
for example in Ref. [3]. The coefficients of the terms are different. Nevertheless, the 
qualitative behaviour of the resulting current might be expected to be similar. 
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9. Appendix 

9.1. Separation of H into one- and two-particle operators 

The total Hamiltonian, Eq. fl2]), with general potential V(r) has the natural additive 
separation 

H = H {int) + H v , 

where H^ nt ^ is the contribution of the electron-electron interaction, Eq. ([3]), and where 
by u Hy" the operator 

i 

is denoted. The separation is obtained from a classification into one-particle operators 
(Hy) and two-particle operators (H^ nt ^), Ref. j3]. I tried to express or verify this 
classification also by the use of the antisymmetric product. 

By methods very similar to the ones applied in the main part of the text one obtains 
for the action of Hy and H^ mt ' on any Slater determinant of one-electron wave functions 

H (int 1 [fa <g> . . . <g> Vjv] = Yl s 9< K fv) [H (int \^ f ® ip g )] ® 

l<f<9<N 

H v [ipi <8 ■ ■ ■ <g> iPn) = sgn(Kf)[H v ipf]® 

l</<iV 

1pK f (2) ® ■ ■ ■ ® lpK f (N)- 

The permutations Kf,Kf g £ Sjy are defined in the same way as in the main text, i.e., 
Kf(l) = f, K/ ff (l) = /, Kf g (2) = g, the rest is irrelevant. 

For the action of the operators on a vector space V(R) generated by an arbitrary 
orthonormal system R of one-electron wave functions one obtains in terms of creation- 
and annihilation operators: 

#MV= ^(H v a)®(c a i;), 

a€R 

Hv^ = ~ \- H(mt) (« ® ^)] ® 

for arbitrary ip e V(R). 

If pv(R) is the orthonormal projection onto V(R), then the operator p V (R)H = 
Pv(R)(H( mt ) + H v ) has the representation 

p V {R)H=- A J2 (7®5|^ (mt) (a®/3)) c\c\c p c a 

a,/3,~/,5£R 

+ (Tf\ H v<x) c\c a . (16) 

a,^£R 
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If the orthonormal system R is chosen as a system of eigenf unctions of the one-electron 
operator Hy, then the second line of Eq. (fT6l) turns into 

V A c f c 

where X a is the eigenvalue of the eigenfunction a of Hy. The latter operator is diagonal 
in the Slater determinants of the one-electron wave functions in R. 

The potential V(r) is caused by the nuclei in the crystal. If the crystal is perfectly 
regular, then the potential V(r) can be written as a sum of a periodic contribution and a 
rather slowly varying non-periodic contribution: For every r G M 3 one can imagine a ball 
with centre r and positive and fixed radius. The contribution of the nuclei within the 
ball is periodic in the lattice. The contribution of the nuclei outside this ball is slowly 
varying and has larger negative values in the centre compared to its values closer the 
fringe. Therefore, the eigenfunctions of Hy might be expected to produce an electron 
density which is not periodic, but rather concentrated in the centre of the crystal. This 
is, at least in a vague sense, inconsistent with the fact that, if the lattice of the nuclei 
is perfectly periodic, then the density of the positive charges is also periodic and not 
concentrated in the centre of the crystal. With this reasoning I would assume that it 
might be an inadequate ansatz to try to find or get close to eigenfunctions of the total 
Hamiltonian by using Slater determinants of eigenfunctions of Hy. Even if these Slater 
determinants would turn out to be good approximate eigenfunctions of H, it would seem 
unlikely that an eigenfunction of the final Hamiltonian (including also the nuclei) can 
be constructed by the use of the Slater determinants. 



9.2. Hybridization 

Alternatively, I wrote for the action of the Hamiltonian, Eq. (j2J), on an iV-electron wave 
function if) with arbitrary N > 2: 

Hip= |u^ + ^[^) + ^) + c(A, + A,)]U 

l<i<j<N ' 1 •? ' ' 

== £ 

l<i<j<N 

With the definition 

H {N) := r^—, + tt~~ — r [V(n) + V(r 2 ) + cA] 

|n — T2\ JS — l 

(an operator acting on two-electron wave functions only) one can represent the action 
of H on an arbitrary Slater determinant of iV one-electron wave functions as 

H [V>i ® • ■ ■ ® Vat] = Yl s 9 n ^fg) [H {N) (^f ® if> g )] ® 

i<f<g<N 
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If R is an orthonormal system of one-electron wave functions, then one obtains for 
i> e V N (R): 

H ^ = \ll [H {N \a <g> (3)} <g> [c^] 

and 

The representation of the action of if in terms of the operators seems 
appealing, compared to the separation into H^ nt ' and f/y. Since any antisymmetric 
TV-electron wave function can be written as an (infinite) linear combination of Slater 
determinants, the description of the action of H on the space of the iV-electron wave 
functions is equivalent to the description of the action of the operator on two- 

electron wave functions. The representations of H obtained in the appendix are 
independent of the shape of the potential V(r). 
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